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ABSTRACT 

Giant gas planets in close proximity to their host stars experience strong irradiation. In extreme cases photoevaporation causes a 
transonic, planetary wind and the persistent mass loss can possibly affect the planetary evolution. We have identified nine hot Jupiter 
systems in the vicinity of the Sun, in which expanded planetary atmospheres should be detectable through Lycr transit spectroscopy 
according to predictions. We use X-ray observations with Chandra and XMM-Newton of seven of these targets to derive the high- 
energy irradiation level of the planetary atmospheres and the resulting mass loss rates. We further derive improved Lya luminosity 
estimates for the host stars including interstellar absorption. According to our estimates WASP-80 b, WASP-77 b, and WASP-43 b 
experience the strongest mass loss rates, exceeding the mass loss rate of HD 209458 b, where an expanded atmosphere has been 
confirmed. Furthermore, seven out of nine targets might be amenable to Lyo- transit spectroscopy. Finally, we check the possibility of 
angular momentum transfer from the hot Jupiters to the host stars in the three binary systems among our sample, but find only weak 
indications for increased stellar rotation periods of WASP-77 and HAT-P-20. 

Key words. X-rays: stars - stars: activity - planets and satellites: atmospheres - planets and satellites: physical evolution - planet-star 
interactions - binaries: general 


1. Introduction 


The discovery of giant gas planets in close proximity to their 
host stars brought the stability of these planets and their atmo¬ 
spheres into question. Orbiting a s close as two stel lar radii above 
the photosphere of the host star (iHebb et al.ll2009ll . the so-called 
hot Jupiters and hot Neptunes are exposed to strong irradiation, 
[n particular, the absorption of extreme ultraviolet (EUV) radia¬ 
tion ionizes hydrogen and heats the atmospheric gas. The result¬ 
ing high temperatures of about 10000 K can support a steady 
expansion of the atmosphere, which manifests itself in a plan¬ 
etary wind. In its formation this wind is not unlike the solar 
wind (lParkeil[l958t IWatson et alJll981h . Smaller planets with 
low densities experience the strongest fractional mass loss and 
eould lose their hydroge n and helium envel opes, evolving to hot 
Super-Earth like planets (ICarter et al.ll2012ll . 

Indeed, expanded atmospheres have been confirmed around 
two hot Jupiters^In_a_stud^_o£^he_s^stem^JID_20945^Vidal- 
Madjar et al. (l2003h discovered a 15% dimming in the line wings 
of the hydrogen Lya line of the host star when transited by the 
hot Ju piter in the system, whereas the optical transit depth is only 
1.5% (iHenrv et al Jl2000l : ICharbonneau et al.ll2000ll . To date, the 
presence of this upper atmosphere has been confirmed by several 
observations m easuring excess absorption in H i, O i, C n. Si m, 
and Mg I lines dVidal-Madiar et al.lb004l:lBallester et^ 2007; 


Ehrenreich et alJ 20081: 


Linskv et al. l2010r Jensen et al. 2012 


Vidal-Madiar et al.l 2013ll . In the second system, HD 189733, 

an expa nded atmosphere was also confirmed by s everal obser¬ 
vations d Lecavelier des Et ang s e t ar]|2010l 120121 : iJensen et alJ 
l2012t[^n-Jaffel & BallesterKOl^fk^nd the system also exhibit s 
an excess transit depth in X-rays (iPoppenhaeger et ^1^1^ . 
Further tentative detections of exc ess absorption in tr ansit obser¬ 
vations of the systems WASP-12 (iFossati et alJl20ldh . 55 Cancri 


dEhrenreich et alJl2012h . and GJ436 (iKulow et al.ll2014ll hint 
that expanded atmospheres could be a common feature in tightly 
bound gas giants. 


Four of these five discoveries succeeded using Lya transit 
spectroscopy. There are two reasons for this: First, the upper at¬ 
mospheres of hot gas giants should consist mostly of hydrogen 
and helium, and second, the Lya line dominates the UV spec¬ 
trum oOowjnass_stM 2 _des£itejnterstellar_absor£tion_(France 


et al. 120131) . Nevertheless, with today’s instrumentation these 
types of investigations of exoplanetary atmospheres are only 
possible in close-by systems with strong Lya emission. 


Although the presence of expanded atmospheres is com¬ 
monly accepted at least in the two standard cases, the origin of 
the Lya absorption signal and the essential mass loss are not 
clear. Both can be affected by at least three processes, i.e., the 
planetary wind, the stellar wind, and the stellar radiation pres¬ 
sure. In general, the spectrally resolved absorption signal re¬ 
veals the fraction of absorbed stellar emission and the absorption 
width. The transit depth can be explained by an opaque planetary 
atmosphere covering a certain fraction of the stellar disk. The ab¬ 
sorption width of about 1 A, corresponding to +100 km s“', can 
either be produced by a neutral hydrogen cloud with the given 
bulk velocity or by a static cloud with sufficient optical depth. 

Koskinen et al. (I2013al f5l support the idea that the escap¬ 
ing atmosphere produces a sufficient neutral hydrogen column 
density to explain the Lya transit observations, although this 
wind proceeds with bulk velocities of only 10 km s'. The transit 
depth of their model atmosphere is consistent with the measured 
signal of HD 209458 b. However, the larger velocity offset of 
-200 km s * in the_absor£tionjignal_ofHD_189733^Lecavelier 
des Etangs et al. l2012h can probably not be e xplained with this 
model. Among others, lEkenback et al.l (l2010l) argue that charge 
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Table 1. System parameters of the complete sample (see Sect.|2j 


Host star Planet 


System 

Sp. type 

Ttfs 

(K) 

V 

(mag) 

(S-F) 

(mag) 

iJ-K) 

(mag) 

d 

(pc) 

(d) 

Age 

(Ga) 

^pl 

(-^jup) 

Mpi 

(Mj„p) 

Ppi 

(g/cm^3) 

Teq 

(K) 

Porb 

(d) 

a 

(AU) 

e 

TD 

(%) 

HAT-P-2 

F8V 

6300 

8.7 

0.46 

0.19 

114± 10 

3.7 ±0.4 

0.4 

1.1 

8.9 

7.3 

1700 

5.6 

0.068 

0.5 

0.5 

WASP-38 

F8V 

6200 

9.5 

0.48 

0.29 

110±20 

7.5 ±1.0 

1.0 

1.1 

2.7 

2.1 

1250 

6.9 

0.076 

0 

0.7 

WASP-77 

G8V 

5500 

10.3 

0.75 

0.37 

93 ±5 

15.4 ±0.4 

1.7 

1.2 

1.8 

1.3 

1650 

1.4 

0.024 

0 

1.7 

WASP-10 

K5V 

4700 

12.7 

1.15 

0.62 

90 ±20 

11.9 ±0.9 

0.6 

1.1 

3.2 

3.1 

950 

3.1 

0.038 

0 

2.5 

HAT-P-20 

K3V 

4600 

11.3 

0.99 

0.67 

70±3 

14.6 ±0.9 

0.8 

0.87 

7.3 

13.8 

950 

2.9 

0.036 

0 

1.7 

WASP-8 

G8V 

5600 

9.8 

0.82 

0.41 

87 ±7 

16.4 ± 1.0 

1.6 

1.0 

2.2 

2.6 

950 

8.2 

0.080 

0.3 

1.3 

WASP-80 

K7-M0V 

4150 

11.7 

0.94 

0.87 

60 ±20 

8.1 ±0.8 

0.2 

0.95 

0.55 

0.73 

800 

3.1 

0.034 

0 

2.9 

WASP-43 

K7V 

4400 

12.4 

1.00 

0.73 

80 ±30 

15.6 ±0.4 

0.8 

0.93 

1.8 

2.9 

1350 

0.8 

0.014 

0 

2.6 

WASP-18 

F6IV-V 

6400 

9.3 

0.44 

0.28 

99± 10 

5.0±1.0 

0.7 

1.3 

10.2 

10.3 

2400 

0.9 

0.020 

0 

0.9 

HD 209458 

GOV 

6065 

7.7 

0.58 

0.28 

50±2 

11.4 

1.5 

1.4 

0.69 

0.34 

1450 

3.5 

0.047 

0 

1.5 

HD 189733 

K0-2V 

5040 

7.6 

0.93 

0.53 

19 

12.0±0.1 

0.7 

1.1 

1.1 

0.96 

1200 

2.2 

0.031 

0 

2.4 

55 Cnc (b) 

KOIV-V 

5200 

6.0 

0.87 

0.58 

12 

42.7 ±2.5 

6.7 

— 

0.80 

— 

700 

14.6 

0.113 

0 

— 

GJ436 

M2.5V 

3350 

10.6 

1.47 

0.83 

10 

56.5 

6.5 

0.38 

0.073 

1.7 

650 

2.6 

0.029 

0.2 

0.7 


Notes. Columns are: name of the system, spectral type, effective temperature, visual magnitude, colors (SIMBAD or according to spectral type), 
distance, rotation period, and gyrochronological age of the host star (see Sect.O; planetary radius, mass, density, equilibrium temperature (cited 
or Teq = Teff sJRst/la), orbital period, semimajor axis, orbit eccentricity, and transit depth (TD). 


RefereIlceSJ_^hedatawerecotnraleduMng_exo£lanetSJorgJWright_et^y|20^J^_anddlejbllowmg^ublications^^^^P;2JB^ros_etalJJ20^ 


Leeuw en ( 120071) ; iP al et alji 20101 ), from v sin i, Te„ varies due to eccentricity (1250 to 2150 K); WASP-38 : [BatTos et al.l l l2011ll : lBrown et al.l 


ll2012il: WASP-77 : lMaxted et alJ ll2oi3^ : W ASP-10:IChristian et al.r(l2009ll:I.Tohrison et alJ ll2009h : ISmith et al.l i2Q09h, Prot from LSP; HAT-P-20: 


iBakosetalJ (IMTh. from LSP; WA SP-8: IOueloz et alj (l201(]h : ICubillos et al.l ll2012h fr om LSP; WASP-80: miaud et alJlIM^. from 


i; WASP-43: iHellieret all lIMTh: WASP - 18: iHellier etal 


ll2009h : iPillitteri el^jlj j201^. P^ot from vsini; HD 209458 Charbonneauet all 


( l2000h : lHenrv et idJd200(]|) : lToiTes et al]i2008h:ISilva-Valiol(120081): HD 189733:lBouchv et and2^:lHenrv & Winnl(l2^:ISouthworthlll201(]|) : 
c:|B 


55 Cnc: lButler et alJ l^l997^ : ^McArthur et alJ ll2004h : lGrav et al.l <2003h : lFischer et alJ ( l2008h : GJ436: lButler et al.l (l2004h :[K 


nutson et al 




exchange between the fast ionized stellar wind and the neutral 
planetary wind creates the hydrogen population causing the Lya 
absorption, but their model requires 40% of Jupiter’s magnetic 
moment to reproduce the redshifted absorption. In principle ra¬ 
dio observations could be used to determine the magnetic field 
strength of hot Jupiters, but most host stars have no t been de¬ 
tec tecynj;adio^)bseryations_soJ^(|siroAia^et^020l3)- In partic¬ 
ular, iLecavelier Des Etangs et af] ( 201 ll) presented observations 
with the Giant Metrewave Radio Telescope resulting in upper 
limits for the meter wavelength radio emission from HD 209458 
and HD 189733. With certain assumptions, the nondetections 
indicate upper limits for the planetary magnetic field strength 
of only few times that of Jupiter. Regarding the field strength 
theoretical estimates are not particularly helpful either, as they 
predict values from insignificant to dominating field strength 
dTrammell et al.l2()Tll) . Furthermore, neutral hydrogen is also ex¬ 
posed to the radiation pressure from stellar Lya emission, which 
can reach several times the gravitational acce leration of the host 
star (iBourrier & Lecavelier des Etangsll2013h . In a collisionless 
regime, and with sufficient neutral hydrogen supply, radiation 
pressure can produce a neutral hydrogen population at velocity 
offsets of -100 kms ^ Larger velocity shifts or redshifted ab¬ 
sorption cannot be explained by radiation pressure. 


Under certain circumstances photoevaporation, charge ex¬ 
change, and radiation pressure can all affect the atmospheres of 
hot gaseous planets and especially mass loss rates. Since each 
mechanism depends on different system parameters, they should 
be distinguishable by comparing the absorption signals from dif¬ 
ferent exoplanets. A common parameter of all three processes 
is the size of the expanded upper atmosphere, which is set by 


the gravitational potential of the planet and by the EUV irradi¬ 
ation. A planet of fixed size with a smaller mass produces not 
only a stronger planetary wind, but also has a larger interac¬ 
tion region for charge exchange and radiation pressure. Besides 
the atmospheric size, charge exchange is mai nly affected by the 
tempe rature and yelocity of the stellar wind (iHolmstrom et al.l 
l2008h . whereas the radiation pressure primarily depends on the 
Lya emission line_strength_oftheJiosLsto_{Bourrier & Lecaye- 
lier des Etangs l2013h . While the two systems with secure de¬ 
tections of expanded atmospheres already show differences in 
their Lyg absorption signal (e.g., comp are IVidal-Madiar et aTI 
l2003l : lLecayelier des Etangs et alJl2012h . a larger sample of de¬ 
tections, increasing the phase space of system parameters, will 
eyentually reyeal correlations of the absorption signal with sys¬ 
tem parameters and help to identify the dominating processes in 
the atmospheres. 


The X-ray luminosity of exoplanet hosts is a crucial param¬ 
eter for both, estimating the mass loss rate of a hot gas giant and 
assessing the detectability of the expanded atmosphere. In the 
energy-limited case, the mass loss rate depends on the radiatiye 
energy input due to hydrogen ionizing emission of the host stars, 
howeyer, EUV emission is mostly extinguished by interstellar 
absorption. Seyeral methods exist to reconstruct this emission 
and p articularly the X-ray lumino sity is a robust and direct proxy 
(e.g.. ISanz-Eorcada et al.ll201lh . because both spectral ranges 
are formed in associated structures in the stellar atmospheres. 
Eurthermore, X-rays constitute a significant fraction of the high- 
energy radiat iye output from actiye host stars like HD 1 89733 
and Corot-2 (IPillitteri et al. 201^ ISchroter et all 1201 ll: Sanz- 
Eorcada et al. 1201 ih . X-ray luminosities are also closely cor- 
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related w ith the Lva luminosity of main-sequence stars ( Lin- 
sky et al. l2013h . Thus, they can be used to predict the signal 
in Lycr transit observations aimed at a detection of expanded at¬ 
mospheres. 

In an effort to increase the number of planets with detectable 
expanded atmospheres, we have identihed a sample of nine hot 
gas giants in the vicinity of the Sun, where predictions yield 
strong Lya emission amenable to transit spectroscopy. Here we 
report on our Chandra and XMM-Newton observations of seven 
targets without prior X-ray observations. The X-ray observations 
are used to determine the total high-energy irradiation of the 
planetary atmospheres. For this purpose, we compare three dif¬ 
ferent methods to reconstruct the EUV emission of host stars. We 
present the first energy-limited mass loss rates based on observa¬ 
tions. The improved Lya luminosities of our targets in combina¬ 
tion with estimates of interstellar absorption reveal the best tar¬ 
gets for future transit spectroscopy campaigns. Finally, we assess 
the possibility of enhanced stellar rotation induced by tidal inter¬ 
actions with the close gas giants in three binary systems among 
our sample. 


2. Target selection 

To identify the most promising targets for Lya transit spec¬ 
troscopy among conhrmed extrasolar pla nets we first selected 
hot, gaseous planets. The recent work by iMarcv et al.l (l2014ll 
has shown that the transition from rocky to gaseous planets with 
large amounts of volatile elements occurs at ~ 2 R(q. Therefore, 
we selected all planets with radii greater than 2 R(q and an orbital 
distance smaller than 0.1 AU to ensure high levels of irradiation. 
Planets with an optical transit depth smaller than 0.5% were also 
excluded to increase the contrast of the expected absorption sig¬ 
nal. 

In a second step, we predicted the strength of the Lya emis¬ 
sion line, using the following relation between effective temper¬ 
ature Teff, stellar rot ation period Prm, an d the Lya line flux in 
ergcm^^ s ' at 1 AU (iLinskv et al.ll2013Tl : 

log(LLy„(l AU))= (1) 

0.37688 -10.0002061 for P,ot = 3 - 10 d, 

< 0.48243 H- 0.0001632 T^g for = 10 - 25 d, 

-1.5963 -H 0.0004732 T^g for > 25 d. 

The stellar rotation periods were derived from the stellar radius 
and rotational velocity (see Sect. 12. lb . Finally, we scaled the line 
flux with the stellar distance to obtain the flux at Earth (see Ta- 
ble|4]i. Interstellar absorption has a decisive effect on the observ¬ 
able strength and profile of the Lya line. Because of considerable 
uncertainty in predicting the absorption, it is neglected in the se¬ 
lection process, but for the final targets we provide estimates in 
Sect. 14.21 The distance scaling ensures that the most promising 
systems are included. 

Based on the anticipated Lya flux, we ranked the preselected 
systems to And the most suitable targets for transit spectroscopy 
campaigns. The ranking order is mostly dominated by the dis¬ 
tance of the host stars so that the most promising targets have 
distances less than 120 pc. Thus, they are probably contained 
in the Local Bubble of low-densitv, hot interstellar gas (Red- 
field & Linsky 120080. This further limits the uncertainty intro¬ 
duced by interstellar absorption. Among our Anal sample we ac¬ 
cepted only targets with a Lya flu x exceeding 1/5 of the unab¬ 
sorbed line flux from HD 209458 (I Wood et al.ll2005l) . Consid¬ 
ering the remaining uncertainty in the prediction process, these 
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Fig. 1. Periodogram of the available stellar photometric data (Super- 
WASP) and the folded light curves, with the adopted stellar rotation 
periods marked by pointers. Periods with a power in excess of 30 occur 
with a rate of less than 1 in 10000 by chance. Substructure of major 
period peaks is caused by period or phase shifts between different ob¬ 
serving seasons. The 1 d period is prominent in all periodograms. 






targets can exhibit sufficiently strong Lya lines for the detection 
of expanded atmospheres. 

With this procedure we found 11 hot gas giants among 
all confirmed planets. As expected, the best-known systems, 
HD 189733 and HD 209458, lead the ranking. The three other 
systems with reported detections of absorption signals are not to 
be found in this target list. The distance of WASP-12 is too large 
for Lya transit spectroscopy; 55 Cncb is not transiting the host 
star, but only the expanded atmosphere has been found to un¬ 
dergo a grazing transit; and GJ 436 b is a slowly rotating host 
star with no certain detection of the equatorial rotational ve- 
locitjyisedto_derive_estimates_Jorflie^^^himinosit^Lanotte 
et al. 20141) . 55Cncb and GJ436b have proven to be suitable 
targets for Lya transit observations, thus, we included them in 
our analysis. The sample was scanned for previous X-ray obser¬ 
vations: lAABPGX_wsno^etected in a dedicated X-ray obser- 
va tion (iPillitteri et al.l l2014l) and WASP-43 has been analyzed 
bv ICzesla et al.l (120131) . The remaining seven host stars did not 
have any prior X-ray observation and unrestrictive ROSAT upper 
limits. The system parameters of our targets are summarized in 
Table □ 


2.1. Stellar rotation periods 

Precise stellar rotation periods of the targets are needed to es¬ 
timate the Lya emission line strength. We compared published 
values with estimates based on the stellar rotational velocity and 
our own periodogram analysis of the host star’s light curves (Su- 
perWASPO). The adopted values are given in Table[T] 

An initial estimate for the period can be derived from the 
stellar radius Rs and the equatorial rotational velocity Vgqi 

Prot = 27rPs X Vgq' , (2) 


' http://WWW.superwasp.org/ 
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Table 2. Details of the observations 


Target 

ObsID 

Inst. 

Start time 
(UT) 

Duration 

(ks) 

Start 

phase 

Stop 

phase 

Transit start 
phase 

HAT-P-2 

15707 

C 

2013-11-1601:15 

9.9 

0.971 

0.996 

0.984 

Wasp-38 

15708 

C 

2014-01-18 07:52 

9.9 

0.984 

1.004 

0.986 

Wasp-77 

15709 

c 

2013-11-09 13:36 

9.9 

0.883 

0.997 

0.967 

Wasp-10 

15710 

c 

2013-11-15 07:31 

9.9 

0.786 

0.830 

0.985 

Hat-P-20 

15711 

c 

2013-11-2417:33 

9.9 

0.699 

0.749 

0.987 

Wasp-8 

15712 

c 

2013-10-23 20:42 

9.9 

0.109 

0.127 

0.989 

Wasp-80 

0744940101 

XMM 

2014-05-15 20:13 

15.5 

0.718 

0.783 

0.986 


Notes. ’C’ is a Chandra ACIS-S observation, XMM is aaXMM-Newton EPIC PN observation. Start and stop phase indicate the exposure duration 
in reference of the planetary orbit with the transit occurring at phase zero. 


assuming that the inclination of the rotation axis is 90° — a rea¬ 
sonable assumption for most transiting systems. The rotational 
velocity and the stellar radius are available for most systems, and 
in the cases of WASP-80, HAT-P-2 and WASP-18 more precise 
rotation periods are not available. 

WASP-38, WASP-77, WASP-10, HD 209458, HD 189733, 
55Cnc, and GJ436 have published rotation periods based on 
photometric variability, which correspond well with the values 
derived from Eg. El The rot ation period of WASP-43 was ana¬ 
lyzed bv lHellier etal.1 (1201 Ih . 

We download the original photometric data of the available 
targets (HAT-P-2, WASP-10, HAT-P-20, WASP-8) from the Su- 
perWASP archive and analyzed th e light curves via a generalized 
Lomb-Scargle periodogram (LSI ^Zechmeister & Kurst^l2009ll 
using the PyAstronomy packagq^ to derive stellar rotation pe¬ 
riods (see FigllJ. WASP-10 and HAT-P-20 show isolated peaks 
with high powers in the periodogram that agree with the v sin i 
based rotation periods within the errors. The analysis of HAT- 
P-2 remained inconclusive, so we reverted to the rotation period 
from Eq. |2 

The data of WASP-8 consist of two seasons with different 
major periods, which have false alarm probabilities of less than 
10 Only the period around 16.4 d is present in both seasons 
and shows a clear photometric variation in the phase folded light 
curve. The estimate from the stellar rotational velocity and the 
radius is about a factor of two larger, but the planet is in a retro¬ 
grade orbit with a projected a ngle of 123° betwee n the orbital 
and the stellar rotation axes (lOueloz et aP 1201011 . Hence, the 
prior assumption of an inclination angle of 90° is probably in¬ 
valid, and the photometric rotation period in combination with 
the V sin i indicate an inclination of the rotation axis with the line 
of sight of about 30°. 

3. Observations and data analysis 

We observed six planet hosts with the Chandra X-ray obser¬ 
vatory using 10 ks long exposures and WASP-80 with XMM- 
Newton for 16 ks (see Table EJ. Three Chandra observations 
partially covered the planetary transit; the exposure of WASP- 
38 occurred completely within the transit of the hot Jupiter. 

We analyzed the Chandra ACIS-S observations with the 
Chandra Interactive Analys is of Observations sof tware package 
(CIAO) 4.6, CALDB 4.5.9 (iFruscione et al.ll2006ll . We extracted 
the source counts in a circular region with a radius of 2" and 
background counts in another circular region with a 25" ra¬ 
dius placed in a source-free region close to the target position. 
The target coordinates at the observational epoch were obtained 

^ http://www.hs.uni-hamburg.de/DE/Ins/Per/Czesla/PyA/PyA/ 


Table 3. Optical Monitor results of WASP-80 


Piker 

Cen. Wavelength 
(nm) 

Plux“ 

Mag. 

U 

344 

5.4 + 0.5 

14.5 

UVWl 

291 

1.5 

16.0 

UVM2 

231 

0.1 

19.2 

UVW2 

212 

<0.04 

<20.3 


Notes. The UVW2 filter did not yield a detection of the target. Count 
rate conversion for a KOV star (XMM-Newton User Handbook). The 
UVW1 and UVM2 conversion introduces a factor of two error and for 
UVW2 a factor of ten. Magnitudes were computed according to the 
fluxes. (10“'^ ergcm“^s“‘A“*) 


from the SIMBAD database, accounting for proper motion. For 
HAT-P-20, we used th e proper motion from the UCAC4 catalog 
(IZacharias et al.ll2013h . 

In Fig.|2]we show the source regions filtered for photon ener¬ 
gies between 0.2 - 2 keV, with the extraction regions marked by 
circles. We detected six host stars as X-ray sources (including 
the XMM-Newton observation, see below), only WASP-38 is a 
nondetection with zero photons at the source position. The cen¬ 
troids of the source counts in the Chandra pointings show offsets 
smaller than 0.5" from the nominal source positions, while the 
Proposers’ Observatory Guide (Version 16) cites a one sigma 
pointing error of 0.6 "q In the WASP-10 exposure, two pho¬ 
tons were detected close to the source position, which shift the 
centroid by 1.1" if attributed to the target, but only 3% of the 
Chandra observations show a pointing error in excess of 1". The 
proper motion correction amounts to only 0.5" with a cited error 
of 10%. Accordingly, it is most likely that these are background 
photons, so they are not included in the analysis of WASP-10. 
In any case, the two photons do not have a large impact on the 
analysis. 

WASP-80 was observed on May 15th, 2014, with all three 
detectors of the European Photon Imaging Camera (EPIC) on 
board of XMM-Newton in full frame mode with medium filters; 
more details about the telescope and the detectors can be found 
in the XMM-Newton User Handboolfl We also obtained four 
full frame images with the Optical Monitor (OM), each 1320 s 
long with the filters U, UVWl, UVM2, and UVW2. The derived 
fluxes and magnitudes are given in Table [3] We reduced the data 
with the Scientific Analysis System (SAS 13.0.0) with the stan¬ 
dard procedure and filters. The X-ray exposures are free of high 
background periods, so that no time filtering was applied. The 

^ http://cxc.harvard.edu/proposer/POG/ 
http://xmm.esac.esa.int 


Article number, page 4 of 1 11 1 
























M. Salz, P. C. Schneider, S. Czesla, J. H. M. M. Schmitt: Irradiation and mass loss of hot Jupiters 


2:28 

38.0 37 

.6 37 

.2 36 

.8 

q 

■ 


■ 

■ 

■ 

:40.0 36 


r? 



o 

q 


V,^ 

X. 









6.8 36 

4 16:2C 

:36.0 35 

.6 35 

.2 34 

.8 

o 

id 

Din 






o ii/ 

T-P-2 





in V 

o 













g 

■ 

.2 58 

.8 58 

.4 23:15 

:58.0 57 

.6 


q 

J/V 

\SP-10 




o 






q 

d 






rsi 

M ■ 






3 6 9 12 15 18 21 24 27 


.0 bo 

4^ 

rj 

.4 7:27 

40.0 3S 

■ 

.6 ■ 3S 

.2 

rH 

q 

ni 

rH 

< 

X 

”x 

■-P-20 



S 

fN 

q 





o 

o 





o 






123456789 


123456789 


36 

q 

c6 

.8 36 

■ 

,4 23:59 

■ 

:36.0 35 

■ 

.6 3! 

'.2 

(—1 
o 

q 


WASP-8 




in 

1 o 






in 

o 

d ■ 

o 


X 




d 







2 4 6 8 10 12 14 16 18 


0.5 1 1.5 2 2.5 3 3.5 4 4.5 


0 44 

.0 42 

'.0 20:12 

■40.0 ^ 

.O'--^ 36 

.0 

■ 

0 ■ 


^ - 

1 

'-■feo “ 

«■ 

■ ■■ 


0 

m ■ 

0 . - B 

d " 
0 

■ ^ ■„ 

^ S - j 


■■ , • 
■/ ■ s 

. ■' 

dO 

(N 

0 



■ ■ 

■ ■■ ■ , J 

■ 


123456789 


Fig. 2. Source regions filtered for soft X-ray emission (0.2 - 2.0 keV). The target positions are marked by circles with a 2"/15" radius in the 
ChandraIXMM-Newton observations; the position of known companions is marked by crosses. In the panel of WASP-80, a nearby detected 
2MASS source is also marked by a cross. The X-ray source close to HAT-P-2 is discussed in Sect. 13.II Six out of seven targets have been detected 
as X-ray sources. 


source counts were extracted in a circular region with a radius of 
15" and background counts were extracted from three circular 
regions on the same CCD, which were chosen with at least 30" 
distance to any recognizable X-ray source. 

3.1. Binaries 

WASP-8, WASP-77, and HAT-P-20 are known visual binary 
stars and the positions of the faint red companions (BD-07 436B, 
2MASS J07273963-H2420171, WASP-8 B) are marked in Fig.|2] 
by crosses. The relative position of the companions was obtained 
from images^f the Two Micron All Skv Survev (2 MAS S. Skrut- 
skie et al. I2006I) and transfered to the X-ray images, assuming a 
common proper motion of the binaries. For the companion of 
HAT-P-20, the infrared colors J = 10.14, H = 9.44, and K = 
9.22 from 2MASS indicate an M-type dwarf. The companions 
of WASP-77 and WASP-8 have spectral types of K 5V and M, 
respectively (iMaxted et al.ll201^IOueloz et al.ll2010h . 

The 2MASS images of all targets were screened for fur¬ 
ther close companions. Only WASP-80 shows a 2MASS source 
(2MASS 20124062-0208333) about 4 magnitudes darker at 9" 
distance, which is also marked in Fig.|2j the source is also clearly 
distinguishable in the OM exposure. The infrared colors indi¬ 
cate a late K to early M-type star, similar to the spectral type of 
WASP-80 (iTriaud et al.l2oi3ll . which is not known to be a binary 
star. The magnitude difference at a similar spectral type suggests 
that the 2MASS detection is a background source. 

In the exposure of HAT-P-2, we detected a soft X-ray source 
at a distance of 16" from the target, corre sponding to 1800 AU 
at the distance of the HAT-P-2 (114 pc. Ivan Leeuwenll2007ll . 
The 2MASS J-band yields no detection with an upper limit of 
15.8 ma^ hence a dwarf with a spectral type later than M5 

^ http://www.ipac.caltech.edu/2mass/overview/about2mass.html 


would not be detected. Assuming the distance of HAT-P-2, the 
X-ray source has a luminosity of 2 x 10^^ ergs"', whereas the 
maximal Lx/Thni = = 10"^ in M-dwarfs results in 4 x 10^^ ergs"* 
(iJames et al.ll2000l) . The light curve excludes a strong flare and 
the gyrochronological age of HAT-P-2 (0.5 Ga in our analysis, 
1.6 Ga in the more detailed analysis of lBrowr3l2014ll excludes 
increased activit^_due_Jo_ajoung_age_(<^100_Ma^Stelzer & 
Neuhauser l200lh . We conclude that this source is most likely 
a background object, possibly of extragalactic origin because 
HAT-P-2 does not lie in the galactic plane. 

3.2. Light curve anaiysis 

For our analysis of the high-energy irradiation of planetary at¬ 
mospheres, we are interested in the mean coronal emission of 
the targets. Strong flares during the observations would result 
in unreasonably high luminosity estimates. The exposure time 
corrected X-ray light curves of the targets with 1 ks binning are 
shown in Fig. [3] We added the background subtracted counts 
of the three EPIC instruments for the light curve of WASP-80. 
No background correction is applied in the Chandra exposures. 
None of the observations shows a strong flare. 

The count rate of WASP-80 declines significantly during the 
first 6 ks of the observation. The bottom panel of Fig.[3]shows the 
hardness ratio (0.65 - 2.0/0.2 - 0.65 keV) of the target. The ratio 
decreases along with the count rate, indicating that this trend 
can be related to prolonged flaring activity at the exposure start. 
However, the count rate around 6 ks is unusually low compared 
to the rest of the observation, so eventually we use the complete 
exposure to derive the mean X-ray flux of the target. 

The observations of HAT-P-2 and WASP-77 partially cover 
the planetary transit. Individual X-ray exposures are not sensi¬ 
tive to the level of absorption expected from expanded atmo- 
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Fig. 3. X-ray light curves of the targets with a 1 ks binning. The XMM- 
Newlon observation of WASP-80 is about 5 ks longer, and the count rate 
is the sum of all EPIC detectors; the error is given by the pipeline. In 
the Chandra observations, error bars represent the 1-sigma Poissonian 
uncertainty on the count rate. Planetary transits during observations are 
marked by shaded areas. The WASP-10 panel shows only the arrival 
times of the five source photons. The last panel depicts the hardness 
ratio (0.65 - 2.0/0.2 - 0.65 keV) in the WASP-80 observation. No pro¬ 
nounced flare occurred during the observations. 


spheres of hot Jupiters (~ 7%. IPoDDenhaegeret al.ll2()T^ . Never¬ 
theless, the periods affected by the planetary transits are shaded 
in Fig. [3] Comparing the count rate during the transit with the 
pretransit rate, we derive an upper limit for the size of an X-ray 
opaque planetary atmosphere assuming a homogeneous X-ray 
surface brightness of the host star. For HAT-P-2 we have 31 in¬ 
transit counts with 32.7 counts expected, which corresponds to 
a maximum radius of an expanded atmosphere of 7.3 Rpi, and 
for WASP-77 5/5.5 in-transit/scaled out of transit counts were 
observed resulting in upper limit of 6.3 Rpi for the planetary at¬ 
mosphere (95% confidence). 

It is unlikely that the nondetection of WASP-38 is caused by 
an expanded planetary atmosphere that covers the complete stel¬ 
lar disk, because the atmosphere would have to be expanded over 
two stellar radii (24 Rpi), whereas usual es timates assume ex¬ 
pansi ons by only a few planetary radii (e.g., IVidal-Madiar et af] 
l2003h . Although, the derived upper limit for the X-ray luminos¬ 


ity of the host star is low, the value is consistent with predictions 
(see TablelUl and with other activity indicators (see Sect. 13.41 1. 


3.3. Spectral analysis and X-ray luminosities 


Spectra were extracted from sou rce and backgr ound regions and 
analyzed with XSPEC V12.7.1 (lArnaudlll99i^ . The three spec¬ 
tra from the individual cameras on board XMM-Newton were fit¬ 
ted in a joint an alysis. We u sed the C-statistics, appropriate for 
low count rates (ICashllT97^ . According to the XSPEC guide¬ 
lines, the data were binned to contain at least one photon per birQ 
with the exception of WASP-10, for which the five source counts 
were fitted using the unbinned data. We used an absorbed, one 
temper ature, optically thin plasma emission model (APEC. Eos - 
ter et al. l2012h with solar abundances (iGrevesse & Sauvalll99^ . 
We added a second temperature component for WASP-80. In this 
case, the spectra contain sufficient source counts because of the 
longer exposure and a five times higher effective area at 1 keV 
of XMM-Newton compared with Chandra. 

Since the absorption component of the model remained ill- 
constrained in the fitting procedure, we fixed the interstellar hy¬ 
drogen column density using the source distance and an aver¬ 
age inteKtell^ji^drogenjdensit^,.o£OJ^,am_^Redfield & Lin- 
sky 2000l) . The resulting spectra and models are shown in Pig.|4] 
where the spectra were binned to a lower resolution for visualiza¬ 
tion. The inserts show the 1,2, and 3 sigma confidence contours 
of the two model parameters: plasma temperature and emission 
measure. These models are used to compute the X-ray fluxes in 


the 0.124 to 2.48 keV (5 to 100 A) spectral range (see Table|4|. 
The detected targets have mean coronal temperatures ranging 
from 2 to 9 X 10^ K and their X-ray luminosities vary by one 
order of magnitude from 0.7 to 8 x 10^^ erg s“'. 

The short X-ray exposures do not yield detections of the faint 
companions in the three binary systems. However, one or two 
photons are detected at the expected positions of the secondaries. 
Assuming similar coronal parameters for both components, we 
derive upper limits for the X-ray luminosity of the B components 
(95% confidence, see Table|5]l. 


3.4. Optical observations 

To check for activity indicators in the optical, we obtained quasi- 
simultaneous optical spectroscopy for WASP-38 and WASP- 
77. In particular, the spectra were taken with the “Heidelberg 
Extended Range Optical Spectrograph” (HEROS), mounted at 
the 1.2 m “Telescopio Internacional de Guanajuato, Robotico- 
Espectroscopico” (T IGRE) at La Luz observatory in Mexico 
dSchmitt et al.ll2014l) . The spectra cover the 350 - 880 nm range 
with only a small gap around 570 nm at a resolution of about 
20000 . 

We observed WASP-38 for 30 min, less than one hour after 
the X-ray exposure, and we observed WASP-77 twice, 5 h before 
and 11 h after the X-ray observation with each exposure lasting 
45 min. Both the WASP-38 spectrum and the combined WASP- 
77 spectrum show a mean signal-to-noise ratio (SNR) of 25 in 
the blue channel, which is suffici ent to derive the Mount Wilson 
S-index (5MWO, IWilsonlll978h . We further converted 5MWO 
into the chromospheric log R jjj. index, which represents the ratio 
of emission in the Ca i i H and K emission lines and th e stellar 
bolometric luminosity (iNoves et al.lll984l : lRuttenlll984ll . 

Eor WASP-38, we obtained values of S mwo = 0.153 + 0.006 
and log Rjjj^ = -4.87 + 0.04, which are consistent with a basal 

^ http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 


Article number, page 6 of 1 11 1 










































































M. Salz, P. C. Schneider, S. Czesla, J. H. M. M. Schmitt: Irradiation and mass loss of hot Jupiters 



0.5 1 1.5 

Energy (keV) 


0.5 1 1.5 

Energy (keV) 


0.5 1 1.5 

Energy (keV) 


Fig. 4. X-ray spectra of the targets. The plasma emission models are depicted by the histograms and the source count rates by diamonds with error 
bars. The;y^ contours for the 1, 2, and 3 cr confidence intervals of the model parameters are shown in the inserts. WASP-80 was fitted with a two 
temperature model; the contours of the four parameters are not shown. With the exception of WASP-10, which only has five source counts, the 
model parameters are well confined. 


chromos pheric flux level an d thus a low level of chromospheric 
activity (iMittag et al.ll2013]l . This finding is also consistent with 
the absence of detectable X-ray emission during our observation. 
For WASP-77, we derived an S-index, 5 mwo^ of 0.338 + 0.011, 
corresponding to a log/?[jj. value of -4.57 + 0.02. This indicates 
a moderate level of chromospheric activity clearly exceeding the 
basal level, which is, again, consistent with the detection of X- 
ray emission. Between the consecutive nights, we found no de¬ 
tectable variability in the Ha line, covered by the red spectral 
channel where the SNR is higher (~ 40). 

4. Planetary irradiation and mass loss 

While EUV emission (100 < A < 912 A) contributes about 
~ 90% of the radiation power of hydrogen ionizing (A < 912 A) 
emission from inactive stars (ISanz-Forcada et ^1201 Ih . above 
400 A it is completely absorbed by interstellar hydrogen. To de¬ 
rive estimates for the planetary mass loss rate, which crucially 
depends on the radiative energy input, we have to assess the EUV 
emission of the host stars. 

4.1. Stellar EUV emission 

Here, we compare three different methods to derive the EUV lu¬ 
minosity of dwarfs either based on the stellar rotation period, the 
X-ray lum inosity, or the Lyg luminosity . The first method was 
introduced iLecavelier des Etangsl (l2007l) . who apply a relation¬ 
ship between the stellar equatorial rotation velocity and the stel¬ 
lar flux in the S2 ba ndpass of the Wid e Eield Camera on board 
the ROSAT satellite dWood et al.lll994t) . The flux in this band is 
further scaled to the full EUV range based on an average solar 
spectrum: 

Few (1 AU) = 4.6 f ergcm-^ s"' . (3) 

\2.0kms / 



Fig. 5. Estimates of ^ell^EUV ernissi on in t h e sarnp le of Pizzolato 
et al. ( 120031) based on stellar rotation (red crosses). X-rays (green solid 
line) and Lya emission (blue dashed line). A linear fit to the rotation 
based estimates is represented by a red dotted line. Circles indicate stars 
with rotation periods shorter than 1 d. 

Second, ISanz-Eorcada et ^ (1201 Ih derive a relation between 
the X-ray (5 to 100 A) and EUV luminosities of main-sequence 
stars, 

logLeuv = (4.80 + 1.99) -H (0.860 + 0.073) log Lx ■ (4) 

Their method is based on a sample of stars with a full recon¬ 
struction of the emission measure distribution from 10"^ to 10^ K. 
The plasma model is then folded with an atomic e mission model 
to derive the spectral energy distribution. Third, iLinskv et al.l 
(12014l) derive the EUV flux in several 100 A wide bands on the 
basis of the Lya flux. Their relations are based on intrinsic Lya 
fluxes, EUVE measurements (< 400 A), and semiempirical mod¬ 
els (400 < T < 912 A. lEontenla et al.ir2014l) . 
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Table 4. Results from the analysis of the X-ray observations and the mass loss analysis for the host stars in the sample. 




X-ray analysis 





Mass-loss analysis 




System 

T 

EM-^ 

Lx'* 

Lx 

j rot 

log Leuv 

Lljya 

Fhya 

p'dhsD 
^ Lya 

Lxuv'' 

logM 

log M™* 

AM 


(10*’K) 



(10^* ergs *) 

(ergs-*) 

(ergs-*) 




(g s-‘) 

(g s^‘) 

(%) 

HAT-P-2 

4 9+1.4 

Qc O+10.3 

5 2+*^-^ 

-^‘^-0.5 

8.2)1;^ 

39. 

28.9 

29.0 

6.5 

0.0- 7.8 

1.3 

9.52 

9.79 

0.004 

WASP-38 

— 

— 

<0.8 

< 1.1 

9.4 

<28.5 

<28.7 

<3.7 

<0.3- 0.4 

<0.29 

<9.42 

9.74 

0.008 

WASP-77 

-5,9 

-1.6 

^“^-0.2 

\ 4+0.3 

^ '^-0.3 

2.1 

28.6 

28.8 

5.6 

1.3- 4.9 

3.3 

10.85 

10.3 

0.29 

WASP-10 

1.8!!:? 

^•^-6.0 


r 9+0.6 
^•^-1.2 

1.8 

28.6 

28.7 

5.7 

0.0- 6.2 

1.2 

9.90 

9.47 

0.01 

HAT-P-20 

c 2+3.6 

4 

^•^-1.5 

1 7+0-2 
-^•'-0.4 

1 0+^-2 
-0.2 

1.1 

28.5 

28.7 

8.9 

1.9-14.6 

1.2 

9.22 

8.75 

0.001 

WASP-8 

2 q+1.4 


o 1 +0.4 
2’-^-0.6 

9 0+0.6 
^•8-0.7 

1.9 

28.7 

28.9 

8.0 

3.5-13.3 

0.46 

9.53 

8.90 

0.01 

WASP-80 

9 9+0.8 
^•^-0.9 

Q Q+1.9 
^•^-1.9 

1 7+2-3 
^•'-0.6 

1 9+0.4 
^•^-0.4 

^‘^-0.2 

o 

1 + 
p p 

1.6 

28.5 

28.7 

10.8 

0.0-33.6 

1.1 

10.54 

10.3 

0.29 

WASP-43 

— 

— 

— 

0.8)“ 

0.76 

28.5 

28.7 

6.2 

0.3- 8.7 

6.6 

10.84 

10.1 

0.21 

WASP-18 

— 

— 

— 

<0.07 

17. 

<28.0 

<28.3 

< 1.8 

<1.5- 2.9 

<0.97 

<9.36 

10.3 

0.001 

HD 209458 

— 

— 

— 

<0.03 

4.0 

<27.8 

28.6 

15. 

22 

<0.11 

<9.94 

10.6 

0.04 

HD 189733 

— 

— 

— 

^•'^-0.4 

2.1 

28.6 

28.4 

60. 

190 

2.1 

10.73 

10.2 

0.28 

55 Cnc (b) 

— 

— 

— 

0.04 

0.31 

27.7 

28.1 

65. 

440 

0.01 

— 

— 

— 

GJ436 

— 

— 

— 

0.01 

0.06 

27.1 

27.7 

35. 

50 

0.06 

8.92 

8.78 

0.18 


Notes. Columns are: name of the system, coronal temperature and emission measure (10^° cm ^); X-ray flux at Earth (0.124-2.48 keV, 
10 ergcm^^ s '); X-ray lu minosity and predict ed luminosity dPizzolato et al.ll200^ : EUV luminosity (100-912 A. ISanz-Eorcada et al.ll2011h : 
reconstructed Lya luminosity dLinskv et al.ll20l^ : reconstructed Lya flux at Earth (10“''* ergcm^^ s“’); minimum and maximum Lya flux at 
Earth after interstell^^so^ionJf^^£l^_^^_erg_cm_^_s_^J^_2^UV_^uxat^l^ietM^_distmc£_^2=_^^_9J^_Aj_J^^_erg_c]gj_^_s_^J^jnas^ 

Eorcada et al. l201 ih and rotation based estimate dEhrenreich & Deserlll201 ih : and the last column is the total fractional mass loss of the planet in 
its lifetime (see text). 

References. X-rav flux: WASP-43 dCzesla et alJl201.'lh: WASP-18 dPillitteri et al.ll2014h: HD 2094.58. HD 18973.3. 5.5 Cnc. and G.T436 t Sanz- 
Forcada et al. |2011|). R econstructed Lyo- flux: HD 2 09458 dWood et al.ll2005h . HD 189733 dBourrier & Lecavelier des Etangsll2013l) , 55 Cnc 
dEhrenreich et alJl2012h . and GJ436 dFrance et al.ll2013h . Measured Lyo- fluxes of HD209458, HD 189733, and GJ436 are mean integrated 
out-of-transit fluxes in the original HST data. 


To establish a comparability of thes e three approaches, w e 
use a sample of main-sequence stars from iPizzolato et ahl (l2003h . 
which all have measured X-ray luminosities, known stellar rota¬ 
tion periods, and determined stellar masses. To apply the relation 
from iLecavelier des Etangsl we derive a stellar equatorial rota¬ 
tion velocity w ith the stand ard mass-radius relation for main- 
sequence stars (iLacvl [197^ and then invert E q. |2] The stellar 
Lya luminosity is needed for the iLinskv et al.l method, but that 
luminosity is not available for the stars of Pizzolato. However, a 
close cotTelation betwee n the X-ray and Ly a luminosities exists 
for main-sequence stars (ILinskv et aLOOlSh . 

log Liya = 19.7 + 0.322 log Lx . (5) 

The equation is a l inear h t to the K5 to F5 stars in the sample 
from ILinskv et al.l (l2013h . We can now use Eq. |5]to derive the 
Lya luminosity, which is in turn utilized to compute the EUV 
flux. This method results in a second X-ray based estimate, and is 
valid for K5 to E5 stars. Among our sample, WASP-80, WASP- 
43, and especially GJ 436 have a later spectral types resulting in 
a larger uncertainty of the derived EUV luminosity. 

Eigure |5] shows the three es timates for the EUV luminosi¬ 
ties of the iPizzolato et al.) (l2003h sample versus the X-ray lumi¬ 
nosities. The rotation based estimates show the expected scat¬ 
ter, which reflects the accuracy of the cotTelation between X-ray 
emission and rotation period. Eor a better visualization, we show 
a linear fit to the rotation based estimates, leaving out stars with 


periods shorter than 1 d because they are in the saturated regime. 
Each method results in significantly different EUV luminosities. 
The best agreement occurs for inactive stars (log Lx = 27.2), 
where the three cotTelations agree within a factor of 2.5 , but for 
highly active stars the estimate of ISanz-Eorcada et al.l is about 
one order of magnitude higher than both other relations. The dif¬ 
ferent results of the three approaches give an idea of the uncer¬ 
tainty in deriving stellar fluxes in this spectral range. The uncer¬ 
tainty is larger than the long-term variability expected for solar 
type stars. Eor example, in the Sun the chro mospheric emissio n 
changes by up to 120% over one solar cycle dWoods et al.l2005l) . 
which is a factor of two smaller than the best agreement between 
the EUV estimates. 

Eor the fo l lowing mass loss analysis we use the relation from 
ILinskv et al.l (l2014l) to derive the EUV luminosity of the host 
stars in our sample (see Table |4]i. The relation is based on a de¬ 
tailed analysis of the EUV fluxes in eight spectral bands, and 
we consider the combination of observations and chromospheric 
models, validated on the basis of solar spectra, to be the best 
currently available option. 

4.2. Improved Lya flux estimates and interstellar absorption 

Using Eq. |5] we derive improved estimates for the Lya lumi¬ 
nosities of our targets. The mean dispersion of this equation is a 
factor of two smaller than our initial estimates for the Lya lumi- 
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nosities based on the spectral type and the stellar rotation period 
(see Sect.|2]). We provide the new estimates in Table|4]along with 
the expected unabsorbed stellar flux at Earth. These values can 
be directly compared with the u nabsorbed flux of H D 209458 
{Fi,ya = 15 X 10 ergcm~^s~*. IWood et alJl2005h . where an 
expanded atmosphere has been unambiguously measured. How¬ 
ever, interstellar absorption reduces the observable emission line 
strength at Earth. The absorption depends on the interstellar neu¬ 
tral hydrogen column density, on the relative velocity of the star 
compared to the moving interstellar clouds, and on the shape of 
the Lya line. 

The structure of the interstel lar medium has been determ ined 
within 15 pc around the Sun (iRedfield & Linsl^ l2008l) . but 
our targets are found at distances >60 pc. Therefore, we use 
Nai measurements of interstellar absorption of stars with simi¬ 
lar lines of sight as our targets to estimate the neutra l hydrogen 
column density. Erom Table 1 of I Welsh et al.l (l2010ll . with Nat 
absorption measurements along 1857 lines of sight, we select 
stars within 8° angular distance of our targets resulting in three 
lines of sight for WASP-80, for example, and up to 16 lines of 
sight for HAT-P-2. The col umn densitie s are scaled linearly from 
the distances of the stars of lWelsh et al.l to the distance of our tar¬ 
gets and then converted into ah^drogen column densitv (Eerlet 
et al. 1985lb ' 

To evaluate the absorption, we constructed a mean Lya 
line profile using a double Gaussian combining a narrow 
(130 km s ') and a broad (400 km s ') emission line compo¬ 
nent, which typically describes the line profile well dFrance et al.l 
l2013h . The line flux is normalized to the values derived from 
Eq.m with the broad component containing 15% of the flux. For 
main-sequence st ars, the Lya line pro files often show a double- 
peaked structure (IWood et alJl2005h . which increases the de¬ 
tectable line flux compared to our estimates. The line width 
of stars also varies, imposing further uncertainty on the esti¬ 
mates (broader lines are less absorbed). The standard Lya ab¬ 
sorption cross section with a Voigt-profile is used for the absorp¬ 
tion with a Doppler-width of 10 km s *. Deuterium absorption is 
included accord ing to a Deuterium fraction in the Local Bubble 
of 1.56 X 10“^ dWood et al.l l2004 |). The relative velocity of the 
absorber is neglected because lWelsh et aT] d2010l) did not publish 
the kinematics of the interstellar medium. This introduces an un¬ 
certainty of a few 10%, with absorbed fluxes up to 40% higher 
for relative velocities of up to +50 km s“'. 

The absorbed line flux is computed for the minimum and 
maximum hydrogen column densities scaled from the surround¬ 
ing measured lines of sight (see Table|4|. With this procedure we 
take into account the uncertainty in the derived column densities 
due to inhomogeneous interstellar absorption, which is about an 
order of magnitude for most targets. For comparison, Table 0] 
also contains the measured line fluxes of the four stars with de¬ 
tected absorbtion signals. The limiting flux for the Space Tele¬ 
scope Imaging Spectrograph on board the Hubble Space Tele¬ 
scope is about 2.5 x 10~'^ e rgcm~^ s“' (upper limit derived for 
GJ 1214. lRunce et al]l2013h . For most targets the estimates re¬ 
main inconclusive. However, the derived line fluxes for WASP- 
38 and WASP-18 are too low for a detection because of the in¬ 
active state of the stars; these systems are thus not well suited 
for Lya transit spectroscopy. Despite the large uncertainty in the 
absorbed Lya flux of WASP-80, this is the closest target among 
the sample and one of the most promising systems for follow-up 
campaigns. 


4.3. Mass loss analysis 


With the computed EUV fluxes, we can now determine the mass 
loss rates from the total radiati ve flux fxuv - Fx + fpriv, which 
impingeson^ie_atmos 2 heresjEAaeve^^ 2007 ^Sanz-Forcada 
et al. l201lh 


M = 


3 riFxuY 

4KGp^x 


(6) 


Here, ppi is the planetary density and G denotes the gravitational 
constant. We adopt a heating efficiency of rj - 0.15, which ac¬ 
counts for the fraction of the radiative energy needed for ion¬ 
ization processes or lost throug h radiative cooling . This c hoice 
agrees with recent results from IShematovich et al.l (l2014ll . who 
find T] < 0.20 in the atmosphere of HD 209458 b. Furthermore, it 
en ables a direct comparison wi th the mass loss rates determined 
by lEhrenreich & Des ^ (1201 i h, who utilized the stellar rotation 
based method from Lecavelier des Etangsl (l2007h to estimate the 
EUV fluxes (see Sect. Si. The parameter /T is a reduction fac- 
tor for the gravitatio nal potential of the planet due to tidal forces 
(lErkaev et al.ll2007h . It is given by 


3 1 

K{f) = 1 -+ —7 , 

^ 2 ^ 2^3 


(7) 


with^ = rRL/rpi « (d/3)'^3^^ where b = Mpi/Mjt is planet to star 
mass ratio, and A - a/r^\ the ratio of the semimajor axis to the 
planetary radius. All results are summarized in Table|4] 

The X-ray based mass loss rates of our targets are on aver¬ 
age higher tlwn the rotation based_gstimates from Ehrenreich & 
Desert, which corresponds to the expected offset due to the dif¬ 
ferent methods of deriving the EUV irradiation. However, based 
on our measured irradiation levels, we find the mass loss rate of 
WASP-38 to be more than a factor of two lower and the value 
of WASP-8 a factor of four higher than the previous estimates. 
The mass loss rate of WASP-18 is even one order of magnitude 
lower than the stellar rotation based estimate. 

We caution that these estimates are based on the assump¬ 
tion of energy-limited atmospheric escape, which in general pro- 
videsanjm£eLJigiMoJhe_^ass^Ioss_^ate_of_^^lane^Watson 
et al. ll981h . Furthermore, for the heating e fficiency, values rang¬ 
ing fr om 0.1 to 1.0 are used in liter ature dEhrenreich & Desert! 
1201 Ih . although the recent study of IShematovich et al. (12014 ) 
for HD 209458 b somewhat restricts this unconfined parameter. 
This adds to a considerable source of uncertainty in the deriva¬ 
tion of the high-energy irradiation (see Sect.|4]i. Hence, the mass 
loss rates must be viewed as order of magnitude estimates, but 
are more reliable than previous estimates without determination 
of the high-energy irradiation level (e.g., lEhrenreich & Des'eri 

IMTb . 


The current fractional mass loss of all planets in the sam¬ 
ple is small. The close proximity and the small planetary masses 
of WASP-77, WASP-43, and WASP-80 result in the strongest 
mass loss rates on a par with the mass loss rate of HD 189733 b. 
WASP-80 loses 0.10% of its mass in 1 Ga, assuming constant 
stellar emission. Moreover, high-energy emission is up to a fac¬ 
tor of 100 stronger in young stars (age <100 Ma) (iRibas et al.l 
I 2 OO 5 I : IStelzer & NeuhauseiiT200lh . Assuming this high irradia¬ 
tion level for the first 100 Ma, the planet would have lost addi¬ 
tional 0.28% of its mass. For our targets, we derive the total mass 
loss estimate by combining this value and applying the present 
mass loss rate for the remaining lifetime (for the age estimates 
see Sect. |5]l. According to these values, six hot Jupiters have lost 
less than 0.01% of their masses through photoevaporation (see 
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Table 5. Age estimates for the binary systems within our sample 


System 


Sp.T 

Sep. 

(") 

log Lx 

(ergs^‘) 

Aiso 

(Ga) 

Agyro 

(Ga) 

Ax 

(Ga) 

^tide/^P 


A 

G8V 


28.1 

5.3 

1.7 

4.5 

0.12 

WASP-77 

B 

K5V 

3.3 

<27.5 



>8.9 


HAT-P-20 

A 

K3V 

6.9 

28.0 

7 

0.8 

1.9 

0.066 


B 

M 


<27.4 



> 1.7 


WA QD Q 

A 

G8V 

A C 

28.4 

<3.6 

1.6 

1.6 

0.002 


B 

M 


<27.6 



>1.5 



Notes. Columns are: system name and component, spectral ty pe, sep¬ 
aration of the component s, X-ray luminosity, isochrone age iBrownI 
(l2014ll : lBakos et aTT(l201lh, gyrochronological age, age estimate based 
on the X-ray luminosity jPoppenhaeger & Wolldl2M4h . and height of 
the tidal bulge in reference to the photospheric scale height. 


Table gli. Even if we conservatively assume an uncertainty as 
large as a factor 100 for the total mass loss estimates, these six 
planets cannot have lost more than 1% of their mass because of 
photoevaporation during their lifetime. 

A hnal remark concerning photoevaporation: HAT-P-2 is the 
strongest X-ray source among our targets, but the mass loss rate 
is small because the planet is massive and has a comparatively 
large semimajor axis. Also the predicted mass loss rates of HAT- 
P-20b and WASP-18 b are very small. These planets exhibit the 
highest densities in our sample. In general the formation of a hy¬ 
drodynamic wind from such compact and massive objects is dis¬ 
putable. Hydrodynamic escape is only possible if the sonic point 
in the escape flow occurs before the exobase level is reached, 
otherwise radiative energy cannot be transformed into a bulk 
flow and only Jean s escape proceed s with a significantly lower 
mass loss rate te.g.. lTian et al.ll2()()5h . The compactness of HAT- 
P-2 b, HAT-P-20b, and WASP-18 b may completely prevent the 
formation of a planetary wind. 


5. Angular momentum transfer from hot Jupiters to 
their host stars 

Giant gas plane ts in tight orbits ra ise substantial tidal bulges on 
their host stars (ICuntz et al.ll200r)t) . If the stellar rotation period 
is longer than the orbital period of the planet, it is possible that 
the tidal interaction induces a tor que on the host star and reduce s 
the stellar spin-down with age dPoppen haeger & Wolk |2014jl. 
Based on measured X-ray luminosities IPoppenhaeger & WolkI 
predict stellar ages of binaries, where the primary features a hot 
Jupiter in a close orbit. Binaries have the advantage of provid¬ 
ing two individual age estimates for the system, and a younger 
age of the primary would indicate continuous angular momen¬ 
tum transfer from the hot Jupiter’s orbital motion. Indeed, the 
two host stars with the st rongest tidal interaction s (HD 189733, 
CoRoT-2) in the study of iPoppenhaeger & Woltd suggest a sig¬ 
nificantly younger age of the hot Jupiter host compared with the 
secondary. In contrast three systems with smaller tidal interac¬ 
tions do not exhibit age differences. 

Th e angular momentum transfer was also studied by iBrownI 
(l2014li . who analyzes the gyrochronological and isochrone ages 
of exoplanet hosts. Reduced spin-down of the host stars due to 
tidal interactions would appear as young gyrochronological ages 
compar ed with the isochrone ages in hot Jupiter bearing systems. 
IBrownI found indications for a general age difference of 1.8 Ga 
between gyrochronological and isochrone age estimates, but no 


correlation with tidal interactions. However, the author does not 
exclude effects in individual systems. We us e the mean of the 
three color based rotation-age relations from IBrownI (1201411 to 
derive gyrochronological age estimates for our targets (see Ta- 

ble[B. 

For the three binaries among our targets we derive X-ray 
based age estimates, using the relationship from Poppenhaeger 
& Wolk (1201 ^* The upper limits for the X-ray luminosity of the 
secondaries provide lower limits for their ages. The strength of 
the tidal interactions can be assessed by the fractional height of 
the tidal bulges in reference to the photospheric pressure scale 
height /itide/JJp (ICuntz et al.ll2000ll . In addition, we provide the 
isochrone ages of the three targets from literature in Table |5] 

For WASP-77 A, we derive an X-ray age of about 5 Ga and 
the B component should be older than 9 Ga. For HAT-P-20 A and 
B as well as for WASP-8 A and B, we derive ages of about 2 Ga 
with only slight differences, but the upper limits indicate that the 
B components seem to be about twice the age of the primaries. 
In contrast to WASP-8, WASP-77, and HAT-P-20 exhibit low 
X-ray luminosities compared with the rotation-based estimates. 
Their gyrochronological age estimates are consistently smaller 
than the X-ray based ages. The hot Jupiter WASP-77 b raises the 
highest relative tidal bulge on its host star, and this is the only 
system with an apparent age difference between the host star 
and the companion. The isochrone age of this system is consis¬ 
tent with the X-ray age of the secondary, which further supports 
an age difference between the two stellar companions. The age 
difference also exceeds the general offset of the isochrone age es¬ 
timates of 1.8 Ga found bv lBrownl The same argument holds for 
HAT-P-20, but here the upper limit for the X-ray luminosity of 
the secondary does not reveal a significant age difference. There 
is no difference in the four age estimates of WASP-8, which ex¬ 
hibits the weakest tidal interactions. At this point the data do 
not provide strong evidence for a transfer of angular momentum 
from the hot Jupiters to the host stars. 


6. Conclusion 

We measured the X-ray luminosities of seven hot Jupiter hosts 
and determined the level of high-energy irradiation and the plan¬ 
etary mass loss rates. Additionally, the two previously analyzed 
targets and four systems with detected atmospheres are included 
in our discussion. According to our estimates, six of the eleven 
planets did not lose more than 1% of their mass as the result 
of a hydrodynamic planetary wind during their lifetime. In our 
sample, WASP-80 b, WASP-77 b, and WASP-43 b experience 
the strongest mass loss rates. Our improved Fya flux estimates 
reveal that in seven of the nine systems expanded atmospheres 
could be detectable through Fya transit spectroscopy by stack¬ 
ing a small number of transit observations. While WASP-80b, 
WASP-77 b, and WASP-43 b are good targets because of their 
strong predicted mass loss rates, WASP-8 b is an interesting case 
because its higher mean density could locate the planet close to 
the transition from a strong photoevaporative wind to a stable 
atmosphere. We And only weak indications for an angular mo¬ 
mentum transfer from the orbiting hot Jupiters to the host stars 
in the two binary systems with expected strong tidal interactions. 

We conclude that currently the systems WASP-80, WASP- 
43, and WASP-77 represent the most promising candidates to 
search for absorption signals of the expanded atmospheres of 
hot Jupiters. 

^ For the two F-type stars we use the original relation from iBamesI 
(l2007h for the age relation based on (B - V). 
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